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Objectives 

The  objective  of  this  project  was  to  develop  an  efficient  and  accurate  Eulerian-Lagrangian  (EL)  method 
based  on  a  Hybrid  high-order  WENO-Spectral  and  WENO-Central  Difference  Navier-Stokes  solver  that 
combines  an  excellent  shock  capturing  and  accurate  solution  of  small  scale  flow  features  and  particle 
physics  with  computational  efficiency.  Another  aim  was  to  validate  the  method  against  published  cases 
of  a  shock-particle  cloud  interaction. 


Approach 


Algorithm  development  concentrates  on  the  high-order  coupling  of  the  Lagrangian  particle  solver  to  the 
Eulerian  Hybrid  WENO-spectral  element  and  WENO-Finite  Difference  scheme  and  the  development  of  a 
parallel,  particle  Navier-Stokes  solver.  The  Hybridization  of  high-order  WENO  methods,  Spectral  and 
Finite  Difference  Methods  ensures  that  shocks  and  small  scales  are  captured  in  a  non-oscillatory  manner 
by  WENO,  while  Spectral  and  Finite  Difference  methods  improve  accuracy  and  computational  efficiency 
in  non-shocked  regions.  The  high-order  and  stable  interpolation  between  the  particles  and  the  flow 
ensures  that  the  favorable  characteristics  of  the  Euler  solver  are  preserved  in  the  high-order  Eulerian- 
Lagrangian  solver.  The  performance  of  the  method  is  assessed  through  published  benchmark  cases. 
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Results 


1.  Hybrid  WENO-Finite  Difference,  Navier-Stokes  Solver 


We  have  developed  and  assessed  the  performance  of  a  parallel  Hybrid  WENO-Finite  Difference  Navier- 
Stokes  solver  [1].  Depending  on  the  ratio  of  WENO  usage  versus  Finite  Difference  usage  in  a 
computational  domain,  the  Hybrid  improved  computational  efficiency  by  sever  factors.  With  an 
increasing  order  of  the  scheme,  the  capturing  of  growth  of  small  scale  instability  improves  (Fig.  1). 
Effectively,  a  lesser  grid  resolution  is  required  with  an  increased  order  of  approximation  that  improves 
computational  efficiency  further  by  several  factors  depending  on  the  problem. 


Figure  l.Left  figure:  Iso-surfaces  of  the  number  in  a  Hybrid  computation  of  a  sonic  air  injection  in  supersonic  cross  stream 
visualizing  both  shocks  and  small  scale  structures.  Right  Figure:  Based  on  a  Multiresolution  analysis,  a  WENO-Z  scheme  is 
used  in  the  red  surface  areas  where  sharp  gradients  occur  including  shocks  and  sharp  shear  layers,  while  a  Finite  Diffrence 
scheme  is  used  in  other  areas  with  small  scale  turbulent  structures. 


The  performance  of  the  Hybrid  was  assessed  for  a  benchmark  sonic  injector  in  supersonic  cross  stream 
(Fig-  2). 


Figure  2.  A  Mach  contour  slice  visualization  along  the  center  surface  of  a  3D  computation  of  sonic  air  injection  in  supersonic 
cross  stream.  A  seventh  order  computation  (right  figure)  captures  instabilities  in  the  shear  layer  emanating  from  the  base  of 
the  sonic  injector  while  a  lower-order  computation  on  the  same  grid  (left  figure)  dissipates  the  instability  entirely. 
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2.  Higher-Order  Eulerian-Lagrangian  Hybrid  WENO-Z-Finite  Difference  based  method 

A  higher-order  Hybrid  WENO-Z-Finite  Difference  based  Eulerian-Lagrangian  method  was 
developed  [2].  A  higher-order  ENO  method  determines  the  stencil  of  interpolation  based  on  a 
smoothnes  determination  trhough  divided  differences  on  the  Eulerian  grid  (Fig.  3).  The  ENO 
interpolation  prevented  spurious  Gibb's  oscillations  that  were  shown  to  yield  inaccurate  and/or 
unstable  results  if  Central  Interpolation  was  employed.  High-order  B-Spline  functions  (Fig.  3) 
weight  the  particle  forcings  to  neighbourings  Eulerian  grid  points  ensuring  smooth  sources  in 
the  Eulerian  conservation  equations  that  account  for  the  coupling  from  the  Lagrangian  to  the 
Eulerian  grid. 


Figure  3.  The  particle  solver  was  coupled  to  the  higher-order  Hybrid  Euler  solver  in  a  higher-order  manner  that  prevent 
spurious  Gibbs  oscillations  through  an  ENO  interpolation  (left  figure)  and  a  higher  B-Spline  function  (right  figure). 


The  performance  of  the  Eulerian-Lagrangian  method  was  assessed  against  a  published 
experiment  of  the  interaction  of  a  cloud  of  bronze  particles  with  a  moving  shock  (Fig.  4)  in  ID 
[2],  2D  [2,3,7]  and  3D  [4],  With  increasing  order  of  the  scheme,  small  structures  were  captured 
more  accurately.  Hybridization  of  Central  and  ENO  interpolation  based  on  the  Multi-Resolution 
coefficient  reduced  computational  cost  by  several  factors. 


Figure  4.  Higher-Order  Eulerian-Lagrangian  computation  of  the  interaction  between  the  accelerated  flow  behind  a  moving 
shock  and  a  cloud  of  bronze  particles  show  that  the  higher  Eulerian-Lagrangina  method  sharply  captures  shock  while 
resolving  small  scale  phenomena  in  the  particle  and  gas  phase. 
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The  code  was  extended  and  validated  to  include  the  two-way  coupled  effects  of  evaporation 
between  a  liquid  particle  and  the  gas  phase  [6]. 

3.  Higher-Order  Eulerian-Lagrangian  Spectral  based  method 

Higher-order  interpolation  methods  were  developed  to  couple  Eulerian  Spectral  methods  with 
Lagrangian  particle  methods  [9,10].  An  interpolation  based  on  the  spectral  basis  ensured  that 
high-order  accuracy  of  the  spectral  solver  is  preserved  in  the  particle  tracer  [9].  Combined  with 
weighting  methods  based  on  B-Spline  projections  onto  a  spectral  grid,  the  Spectral  based 
Eulerian-Lagrangian  algorithm  spectrally  converges  to  analytical  solutions.  A  weighting  function 
was  developed  that  is  based  on  a  mixed  polynomial  approximation  to  a  delta  function  projected 
onto  a  spectral  basis  [10].  This  weighting  function  ensured  an  accurate  capturing  of  singular 


interfaces  between  particle  clouds  and  gasses,  while  providing  higher-order  convergence  to 
exact  solutions  away  from  the  interface  (Fig.  5). 


Figure  5.  The  sharp  interface  induced  in  an  Eulerian  solver  by  a  singular  particle  (left  figure)  is  captured  accurately  if  a  the 
weighting  function  is  approximated  by  a  mixed  polynomial  representation  of  a  delta  function,  while  higher-order 
convergence  is  achieved  away  from  the  interface  with  increasing  degrees  of  freedom,  N  (right  figure). 


4.  Parallelization 

The  Eulerian  and  Eulerian-Lagrangian  solvers  were  both  parallelized  based  on  spatial  domain 
decomposition  and  partitioning  with  an  acceptable  scale-up.  Expensive  subroutines  in  the 
Hybrid  solver  were  ported  to  GPU  and  show  moderate  acceleration.  A  high  degree  of 
fragmentation  of  the  data  with  WENO  coverage  for  certain  type  of  shocked  flows  (Fig.  1), 
reduced  the  vectorization  efficiency  and  increased  cache  misses  and  the  frequency  of  memory 
dumps.  An  active  area  of  research  is  to  scale-up  the  Hybrid. 
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